Black hole quantum tunnelling and black hole entropy correction 
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Parikh-Wilczek tunnelling framework, which treats Hawking radiation as a tunnelling process, is 
investigated again. As the first order correction, the log-corrected entropy-area relation naturally 
emerges in the tunnelling picture if we consider the emission of a spherical shell. The second 
order correction of the emission rate for the Schwarzschild black hole is calculated too. In this 
level, the result is still in agreement with the unitary theory, however, the entropy of the black 
hole will contain three parts: the usual Bekenstein-Hawking entropy, the logarithmic term and 
the inverse area term. In our results the coefficient of the logarithmic term is — 1. Apart from 
a coefficient, Our correction to the black hole entropy is consistent with that of loop quantum gravity. 
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In 2000, Parikh and Wilczek proposed an approach to calculate the emission rate at which particles tunnel across 
the event horizon[l|. They treat Hawking radiation as a tunnelling process, and the WKB method is usedQ, In 
this way a corrected spectrum, which is accurate to the first order approximation, is given. Their result is considered 
to be in agreement with an underlying unitary theory. Following this method, a lot of static or stationary rotating 
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black holes are studied Q, Q, fl fl Q fl Q, Q Q Q, Q, H H Q, £ 

The same result, that is, Hawking radiation is no longer pure thermal, unitary theory is satisfied and information 
is conserved, is obtained. But in all of these literature, the entropy of the black hole only contains the Bekenstein- 
Hawking entropy. Will the emission process still consist with the unitary theory if the quantum correction of the 
entropy is taken into account? In present, about the quantum correction, there are different results corresponding to 



different models and methods 
entropy is 4J, |45 1 
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43j . The general formulation of the black hole 
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where a is a model-dependent (dimensionless) parameter. In the case of Loop Quantum Gravity a is a negative 
coefficient whose exact value was once an object of debate (see e.g. but has since been rigorously fixed at 

a = —1/2. In String Theory the sign of a depends on the number of field species appearing in the low energy 
approximation 36]. Therefore, a very interesting work is to introduce the log-corrected entropy-area relation in the 
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tunnelling framework. Moreover, if the emission rate is calculated to the second order approximation, will the entropy 
contain the inverse area term as given in equation ([1])? In this paper, we first show that in the tunnelling picture, 
a logarithm correction term occurs in the expression of the black hole entropy if we take the emission particle as 
a spherical surface wave (spherical shell). Therefore, we expect that Parikh-Wilczek tunnelling framework, in fact, 
is in agreement with the unitary theory to the first order log-correction of the black hole entropy. Then, we verify 
that, if we calculate the emission rate to the second order approximation with Parikh-Wilczek tunnelling framework 
and suppose the emission process to be still in agreement with the unitary theory, the entropy of the black hole will 
contain three parts: the usual Bekenstein-Hawking entropy, the logarithmic term and the inverse area term. Finally, 
we give two comments to the Parikh-Wilczek framework and our calculation. 



II. BLACK HOLE TUNNELLING AND THE FIRST ORDER CORRECTION TO THE BLACK HOLE 

ENTROPY 

As mentioned above, Parikh and Wilczek applied the WKB approximation to calculate the emission rate of a 
tunnelling particle (S-shell). We start with a brief review of the WKB method and the barrier penetration. For 
a massless particle (massless shell), because of the infinite blueshift near the horizon, the characteristic wavelength 
of any wavepacket of the S-wave (see is always arbitrarily small there, so that the geometrical optics limit 

becomes an especially reliable approximation. The geometrical optics limit allows us to obtain rigorous results directly 
in the language of particles. That is, the WKB method and the expression of the emission rate are the same as that of 
a classical massive particle. So, in the following discussion we only study the tunneling process of a massive particle 
(massive shell). 

Schrodinger's equation for the motion of a particle in a centrally symmetric field is 

A^jj + (2m/(h) 2 ) (E-U(r))i> = 0. (2) 
Let us consider the following radial equation: 



By the substitution 



1 d , 1(1 + 1) „ 2m . „ „, „ 
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R(r) = X(r)/r (4) 



equation is brought to the form 

d 2 X 2m 1(1 + 1) 



dr 2 +[^ E ~U(r))-^- r ^]X = Q. (5) 
For S-wave, / = 0, the equation of X(r) is 

^ + 2 -^(E-U(r))X = 0. (6) 

Note that, in the Parikh-Wilczek framework, to calculate the self- gravitation reliably the tunnelling particle is consid- 
ered as a spherical shell (S-wave). In this way, when it emits from the black hole the matter-gravity system transits 
from one spherical state to another. So, the de-Broglie wave function of the emission spherical shell should be 



1&(r) = X(r)/r. (7) 
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That is, the WKB wave function of a particle can be written as 



<Kr)=X(r)/r = iexp[^], (8) 

where 

S(r) = S Q (r) + (?)Si(r) + (-) 2 S 2 (r) + ■■■. (9) 
i i 

Substituting (JSJ) into Schrodinger Equation © yields 

S ° = ± j Pr dr - ( 10 ) 

26^ + SS = 0, (11) 

2S' S' 2 + (S[) 2 + S" = 0, (12) 

where we use a prime to denote differentiation with respect to r. 

To evaluate the probability of a particle passing through the barrier, we divide the whole region of motion of the 
particle by two tunnelling points a and b into three parts: ingoing and reflecting region I, barrier region II and the 
outgoing region 1H. The particle moves as a free particle in region I and M, but region II is classically inaccessible. 

In region I, we take the WKB wave function as follows [46[ 



2 1 f a 7T 

Xl (r) = ^am[ H jf Prdr + -] 



i ■ pa ■ ■ pa 

1 , r Z / ITT I I ITT 

= H^ {6Xp[ ft X ^dr + T ]-exp[-- Jf p r dr--}}, (13) 

where w is the velocity of the tunnelling particle. In region E, the WKB wave function is a linear combination of real 
exponentials. Considering the connexion between the oscillating and the exponential solutions at r = a, the WKB 
wave function in region II can be written as 



i , i , r b , „ , i 



X u {r) = — exp [-- | / p r dr |]exp[-- | / p r dr\\. (14) 



And the WKB wave function in region II is 



The probability of barrier penetration is 



X hi (t) = — 3_exp[-i| / p r dr |] exp P / p r dr + : - 1 ]. (15) 



r _ jmu _ v\ip ou t\ 2 _ v(X out (b)/b)' 2 _ 2ImS 

P ~ Jin _ t# m | 2 ~ «(^i„(0)/0)a " 62 '''^l ' 



Let's now calculate the phase space factor corresponding to the black hole tunnelling. For Schwarzschild black hole, 
the line element in Painleve coordinates is 



ds 2 = -c 2 (l - ^)dt 2 + 2c\/^dtdr + dr 2 + r 2 (d0 2 + sin 2 Odcf), (17) 
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and the radial null geodesies are 



dr , 2MG . , s 

dt= ±c(i -v^- } - (18) 



with the upper(lower) sign in Eq. (I18p corresponding to outgoing(ingoing) geodesies, under the implicit assumption 

that t increases towards the future [4 71]. 

But in this paper we consider the tunneling of the massive particle. That is, the outgoing particle is a massive shell 

(de Broglie s-wave) . The massive quanta doesn't follow radial-lightlike geodesies (IT8l) . Similar to Ref. [19( , we treat 

the massive particle as a de Broglie wave and obtain the expression of r. Namely, 

1 l ff00 1 cV - 2M.Gr 

r = v n = —v„ = = , — . 19 

P 2 9 2 .gel 2r V2M& V ' 

Note that to calculate the emission rate correctly, we should take into account the self-gravitation of the tunnelling 

particle with energy oj. That is, we should replace M with M — u> in (jTTJ) and (TiU)) to describe the motion of the 



particle correctly 



BflB. 



The canonical momentum p r and the imaginary part of the action ImSo can be easily obtained. Namely, 

f Pr , . f dH h 

Pr = J dp r = J — = r ' ' 20 ^ 

I mSo =£p r dr=-^-||]. (21) 

The probability of barrier penetration is 

rf , 2ImS 0l r ,Af At , , As As . . 

r P = ^ ex P [-^] = ex P - m ^ - (i | - m ^)], (22) 

where ^ = In this paper, we investigate the transition of the matter-gravity system from one spherical state to 
another at the same energy. This transition corresponds to the production and barrier penetration of the massive 
spherical shell (or massless shell). That is, this process contains two stages. The first stage is the production of the 
spherical shell from the vacuum fluctuation near the event horizon. The second stage is the barrier penetration. The 
rate of transition from the initial spherical state to the final spherical state is therefore 

r(i - /) = r, ■ r p = r„ ■ exp - in - (^ - in ^)]. (23) 

Let's Compare (f23")l with the unitary result in Quantum Mechanics, T(i — > /) =| Mfi | 2 -(phase space factor), which 
is given in Ref. [l|. | M/, | 2 is the probability amplitude of the process, in this case it is related to the production 
rate of the particle in the vacuum fluctuation near the event horizon. Thus, we obtain 

phase space factor = exp [(^ - In ^) - (^ - In ^)]. (24) 



If we bear in mind that 



N f e f 9 -<? i 
phase space factor = — - = = e f % (25) 

Ns e bi 



we naturally get the expression of the black hole entropy to the first order correction 
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III. SECOND ORDER CORRECTION TO THE BLACK HOLE ENTROPY 

Let's now calculate the tunnelling rate to the second order approximation. In order to get the second order correction 
of the black hole entropy, we write the WKB wave function to the second order approximation. Namely, 

X(r)=exp[^M2 + Si(r) + ?S 2 (r)], (27 ) 



where 



^ i r JS % S ' ) dr - (28) 



Like the treatment in section E, the wave function in region I can be taken as 
Xj(r) = ^si n[ i(/% r dr-^ 2 (r)) + 5 

= ^{exp[^(jy r dr~h 2 S 2 (r)) + ^}-e X p[-^(£p r dr-h 2 S 2 (r))-^}}. (29) 
In this region the expression of S 2 (r) is 

In order to reduce to the first order approximation case, the connexion between the oscillating and the exponential 
solutions at r — a should be 



(31) 



sin [\{J\ dr - h 2 S 2 (r)) + 5^4 exp jf \ p r \ dr - h 2 S 2 {r))]. 

r < a r > a 

On the right hand of the connexion (|3Tj) , the expression of S 2 (r) is 

S2 = £ JS \ + S , S " [ l dr. (32) 



The connexion at r = b is 



^exp[-(| / p r dr|-n a 5 2 )]^~=e Xp [i(/ Pr dr - h 2 S 2 ) + -], (33) 



r < b r > b 

and the wave function in region IH is 



11 i r iir 

X m {r) = -^exp[--(ImS -? l 2 ImS 2 )]exp[-(/ Pr dr - h 2 S 2 ) + -], (34) 

y/v h h J b 4 



where 



ImS 2 =Im [ b - {S \t S " ] dr. (35) 

Since 

ip{r) = X(r)/r, (36) 
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in region I, the ingoing flux density is 



Jin = '^"(V'in^^n ~ = v \^in\ = ^2> ( 37 ) 

and in region M the outgoing flux density is 

iout = ~2^^ out ^™t ~ i>lut^out) = v\ip„ ut \ = ^ exp [-|(ImS - fi 2 ImS 2 )]. (38) 

Therefore, 

a 2 2 

T p = jout/jin = ^2 exp [--(ImS - ft 2 ImS 2 )]. (39) 
For Schwarzschild black hole tunnelling, in classically inaccessible region, we have 

S' =p r = -in^r, Sq = -Mr-j, (40) 



and 



From (|4Tj) we can easily obtain 



So, 



"p v 



s ! = "lf -S- sf = i' (41) 



*— 4(* +s n- ( sS»-?- (42) 



77 o- /2 /2 

5 2 dr = ^(^-^). (43) 



Ah A f A 

Substituting (|2"T1l . (j43 f into ([39| and considering T(i — > /) =| M/j | 2 -(phase space factor), yields 



phase space factor = exp [(^ - In ^ + §i) - ( * - In * + (44) 



Comparing ([44]) with (|25| , we get the expression of the black hole entropy to the second order correction 

s « = t-'"t + ^ +c °"*'-' (45) 

which is consistent with an unitary theory and is also in agreement with the general formulation of the black hole 
entropy. The emission rate is 

r(z-/)~e AS <>. (46) 



IV. CONCLUSION AND COMMENTS 



We showed how a log-corrected entropy-area relation can emerge in the tunnelling picture if we consider the emission 
particle as a spherical shell. We also showed that, if the emission rate is calculated to the second order approximation, 
the black hole entropy will contain three parts: the usual Bekenstein-Hawking entropy, the logarithmic term and the 
inverse area term. In our calculation the logarithmic term and the inverse area term is the consequence of requesting 
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the process to satisfying the unitary theory. Apart from a coefficient, Our correction to the black hole entropy is 
consistent with that of loop quantum gravity. In the following, we give two comments to the Parikh-Wilcek method 
and our calculation. 

1) In this paper, we only take into account the emission of the massive particle. The motion of a massless particle 
(S-wave) is very different from that of a massive particle. However, as mentioned in the first paragraph of the section 
II, a massless shell can be treated in the language of particle. That is, for massless shell we can also apply the WKB 
method and obtain the same functional form of emission rate as that of the massive particle. So, Eqs. (|45p and (|46|) 
are also suitable for massless particle's emission. 

2) In the first order approximation, the previous expression of the emission rate can be written in the following 
explicit form 

T ~ exp (AS q ) = (1 - exp (-8nGMcu(l - ^)). (47) 

In Refs. Q and JjJ the authors pointed out that the coefficient of the log-corrected term in the black hole entropy 
should be positive, otherwise, the probability of emission will diverge when the emission particle's mass, to, approaches 
to M. In fact, if we consider the applying condition of the WKB method, the emission particle's mass, to, will never 
approach to M, it is far smaller than the black hole mass M. Here is our derivation. 
The WKB method is established in the conditions: 

h\S;\«\S' a % (48) 

and 

2h\S' S' 1 \«\S 2 \. (49) 
From (|4TJ)) and (jlTj) , above conditions (jig)) and (|4*9"|) can be incorporated into an inequality, that is 

« 1?% (50) 

Considering p r — —inr and 2(M — uu) ^ r < 2M, (f50|) becomes 

2(M-w)»4/-. (51) 

V 7T 

That is, 

M > lu. (52) 

It means that the Parikh-Wilczek framework is only suitable for the emission of the particle whose energy is far less 
than the mass of the black hole. Most of the time in the evaporation of the black hole this condition is satisfied, 
that is, the mass of the emission particle will not approach to the black hole mass M. Therefore, the coefficient of the 
log-corrected term in the black hole entropy is not constrained to be positive. However, in the last stage of evaporation 
the emission will be very strong, and the mass of the emission particle will be very great, the WKB conditions will 
not be satisfied, then one would have to resort to other mechanisms to describe the last stage of the evaporation. 
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